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Abstract. Cantor's famous construction of the real continuum 
in terms of Cauchy sequences of rationals proceeds by imposing a 
suitable equivalence relation. More generally, the completion of a 
metric space starts from an analogous equivalence relation among 
sequences of points of the space. Can Cantor's relation among 
Cauchy sequences of reals be refined so as to produce a Cauchy 
complete and infinitesimal-enriched continuum? We present two 
possibilities: one leads to invertible infinitesimals and the hyper- 
reals; the other to nilpotent infinitesimals (e.g. h ^ infinitesimal 
such that — 0) and Fermat reals. One of our themes is the 
trade-off between formal power and intuition. 
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1. Introduction 

In a recent issue of The American Mathematical Monthly, K. Hrbacek 
et al have argued that analysis needs better axiomatics than Zermelo- 
Fraenkel set theory with the axiom of choice (ZFC), see [T7]. They 
propose a new axiomatic framework that naturally includes numbers 
they call ultrasmall (i.e., infinitesimal). They mention [T7t p. 803] 
that these ideas are not entirely new, and provide a list of references 
the earliest of which is E. Nelson's 1977 text I32], where the author 
outlined his enrichment of ZFC known as Internal Set Theory (1ST). 

While the axiomatic approach has much to recommend itself, we feel 
that a prerequisite for new axiomatics is a good understanding of the 
mathematical structure that stands to be axiomatized, and not vice 
versa. To make a convincing case in favor of a new axiom system, one 
first needs to explain the basics. In the case of infinitesimal-enriched 
continua, the basics amount to understanding the ultrapower construc- 
tion. One of our goals in this text is to give an accessible explanation 
of the latter in the context of Cauchy sequences, as well as providing 
possible alternatives. 

Cantor's completion of the rationals resulting in the field of real 
numbers proceeds by quotienting the space Cq C of all Cauchy 
sequences of rational numbers by Cauchy's equivalence relation. Sim- 
ilarly, the collection C C of all Cauchy sequences of real numbers 
projects to the Archimedean continuum M: 

C^M. (1.1) 
The corresponding equivalence relation ~c defined by 

M ~c w if and only if lim \un — Vn\ =0, 

n— j-oo 

"collapses" all null sequences to a single point € M. Is there another 
way to define an equivalence relation ~ on the space C that would allow 
some null sequences to retain their distinct identity? In other words, 
can one refine Cantor's equivalence relation among Cauchy sequences, 
in such a way as not to "collapse" all null sequences to zero? The idea 
would be that, relative to a new equivalence relation ~, a null sequence 
of reals would become an actual infinitesimal. In other words, we are 
searching for a new notion of "completion", with respect to which the 
real field M can be completed by the addition of infinitesimals. What 
one seeks is an intermediate stage, *Mf := C/ ~, in the projection (II. ip . 
The subscript "f" in the symbol *Mf stands for "finite" (i.e., there are 
no infinite numbers). Such an intermediate stage *Mf would represent 
an infinitesimal-enriched continuum as in Figure 11.11 
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Figure 1.1. Factoring Cantor's map C — t- M through an 
intermediate stage *]Rf. 

Here, if [u]^ is the new equivalence class of a sequence u, then the 
function 

St : *Mf M, 

defined by 

st([M]^) := lim Un G M 

n— >+oo 

is the usual hmit of a Cauchy sequence u = G C. This function 

represents the standard part of [u]^ G *]Rf, that is a standard real num- 
ber infinitely close to the new number [m]^ G *]Rf. The most natural 
way to obtain a ring structure on *]Rf is to define the equivalence rela- 
tion ~ so that it preserves pointwise sums and products. Therefore, we 
expect *Rf to be a ring rather than a field, because it cannot contain 

the pointwise inverse ( — ) of an infinitesimal \un\, since the inverse 

is not a Cauchy sequence. 

In this text, we will explore two possible implementations of these 
ideas. 

2. A POSSIBLE APPROACH WITH INVERTIBLE INFINITESIMALS 

To implement the ideas outlined in Section [H a possible approach 
is to declare two Cauchy sequences u, v E C to be equivalent if they 
coincide on a "dominant" set of indices in N: 

u V <^=^ {n G N I u„ = Vn} is dominant. (2-1) 

For simplicity, we will use the symbol [u] for the equivalence class [n]^ 
generated by u G C 

What is "dominant"? A finite set in N is never dominant; every 
cofinite set (i.e., set with finite complement) is necessarily dominant, 
and we also expect the property that the superset of a dominant set 
is dominant, as well. Moreover, we expect the relation (12. ip to yield 
an equivalence relation. In particular, the validity of the transitive 
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property for generic Cauchy sequences implies that the intersection of 
two dominant sets is dominant. In fact, let us assume that 

\/u,v,wEC: M~f A v^w =^ u^w. (2.2) 

Then, if sets A and B of indices are dominant, it suffices to takqj 



_1 line A jl line B 

"""■■"^1-;^ ifneN\A + iinen\B 

to obtain n ~ 1 and 1 ~ ty, so that u ^ w from (12. 2p . It follows that 
the set {n G N I M„ = Wn} = A (1 B is dominant. Conversely, if our 
family of dominant sets is closed with respect to finite intersections, 
then the relation ~ is an equivalence relation. For example, the family 
of all cofinite sets 

jr:= {5 C N I N\5 is finite}, 

the so-called Frechet filter, satisfies all the conditions we have imposed 
so far on dominant sets. These conditions define the notion of a filter 
on the set N (extending the Frechet filter). 

It is easy to prove that the equivalence relation ~ preserves pointwise 
operations 

[u] + [v] := [{Un + Vn)^^^] and [u] ■ [v] := [{Un ■ Vn)^^^] (2.3) 
so that the quotient 

*Rf := C/~ (2.4) 

becomes a ring. Moreover, the real numbers M are embedded in *Mf as 
constant sequences. Whether or not *Mf is an integral domain depends 
on the choice of the filter of dominant sets. Thus, the product of 
sequences u and w given by 

{0 if n is even I ^— if n is even 

-, Wn '■= S 

if n is odd 10 if n is odd 

is zero, but whether or not [u] is zero depends on whether the set of 
even numbers is considered to be dominant or not. We will solve this 
problem later. 

To show that the relation as in (12. ip is a refinement of the usual 
Cauchy relation assume that u, v & C coincide on a dominant 
set A. Then we have u„„ — v„„ = for some subsequence cr : N — t- N 
(enumerating the members of the set A) . It follows that u v since u 
and V converge. Of course, the relation ~ is a strict refinement because 



"'^Let us note that in N we have e N. 
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if we take Un = and Vn = 0, then u but {ra G N | u.„ = Vn} = 
is the empty set, which is never dominant. 

Whether the idea expressed by the notion of a dominant set as 
in (12. ip can be considered "natural" or not is a matter of opinion. 
An alternative approach would be to define a new equivalence relation 
in terms of the rate of convergence of the difference u — v. A thread 
going in this direction will be presented in Section [6l but here we will 
continue with the approach based on ( 12. ip . If one accepts this idea, 
then it is also natural to define an order, by setting 



yields a further condition on dominant sets. In fact, if A is dominant, 
then defining 



we have that A is dominant if the first alternative of (12. 6p holds; oth- 
erwise N \ A is dominant. A filter satisfying this additional condition 
is called a free ultrafilter^ Using this additional condition, we are now 
also able to prove that *Mf is an integral domain. 

Theorem 1. is an integral domain. 

Proof. Given nonzero classes [u] ^ and [v] ^ 0, both of the sets 
{ra G N I n„ 7^ 0} and {n E N \ Vn 0} are dominant. Therefore so is 
their intersection. □ 

Given an integral domain, we can consider the corresponding field of 
fractions *]Rfrac- Since *]Rf is also an ordered ring, the order structure 
extends to the quotient field of fractions in the usual way. 

Remark 2. In a classical approach to nonstandard analysis, the equal- 
ity on a dominant set (Formula (12.11) ) is applied to arbitrary sequences, 
rather than merely Cauchy sequences. Nonetheless, our field of frac- 
tions is isomorphic to the full hyperreal fielci *M of nonstandard 

^For the sake of completeness, we say that U (=V{I) is an ultrafilter on the set 
I \i% ^U\U IS closed with respect to finite intersections and supersets; and X eU 
or I\X gU for every subset X of /. 

■^Namely, the field obtained as the quotient of using the same ultrafilter; in 
general, the result will depend on the ultrafilter. 




(2.6) 



(2.5) 
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analysis through 

'^frac ^ ( I £ 

where [{qn)n]u is the equivalence class modulo the ultrafilter U. To 
prove this, note that every sequence q G can be written as g = ^ 
for two null sequences u, v. Thus, we can set := g„ ■ (y^+i) 

3. A FREE ULTRAFILTER, ANYONE? 

Our intuition yearns for meaningful examples of free ultrafilters. 
Such can be obtained by using Zorn's lemma. It is possible to prove that 
some weaker form of the axiom of choice is necessary to prove the exis- 
tence of a free ultrafilter. The use of this axiom in modern mathematics 
is routine. Thus, one of its standard consequences is the Hahn-Banach 
theorem, of fundamental importance in functional analysis]^ Yet, one 
consequence of exploiting this axiom is that we don't possess detailed 
information about how free ultrafilters are made. Moreover, this also 
implies that it is not so easy to prove the existence of a free ultrafilter 
satisfying some given and potentially useful conditions. 

Theorem 3. The Frechet filter can be extended to a free ultrafilter. 

Proof. See Tarski [M] (1930). □ 

We have to admit that it is not so easy to judge the idea represented 
by formula (12. ip . Indeed, starting from our definition of pointwise 
operations and order, one can easily guess that this idea is formally 
very powerful. For example, it is almost trivial to extend to *Mf the 
validity of general laws about real numbers, such as the following law: 

Vx, ?/ G M : sin(x + y) = sin(x) cos(?/) + cos(x) sm{y). (3.1) 

In fact, we can extend trigonometric functions pointwise. Namely, we 
extend sin : M — M to *sin : *Rf — > *Mf by setting 

*sin([M]) := [(sin(M„))„gj^] . 

Finally, the law (13. ip extends to *sin and *cos because the set of indices 
n G N where it is true is all of N. We will deal with this extension of 
laws from M to *]Rf in more general terms in Section [5l 

On the other hand, whatever will be the example of ultrafilter we 
will be able to present, it doesn't seems sufficiently meaningful why the 



^Luxemburg |29j explored possibilities of constructing extensions that rely on 
Hahn-Banach only (rather than full choice). 
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infinitesimal 



(-1)' 



should be considered positive rather than 



n.+l 

negative, or vice versa i 

To summarize, the idea of requiring sequences to coincide on domi- 
nant sets, even if it may seem initially forbidding from an intuitive point 
of view, appears to be formally extremely powerful. As an alternative, 
in Section [6] we will present another idea, which is intuitively clear 
but which doesn't seem equally powerful. Which thread one wishes to 
follow would depend on applications envisioned. 

4. Actual infinitesimals, null sequences, and standard 

PART 

What are, formally, the infinitesimals in the ring *]Rf of (12. 4p . and 
how are they related to null sequences? An infinitesimal is a number 
in *Rf which belongs to every interval of the form : 

Definition 4. We say that x G *]Rf is infinitesimal if and only if 

Vn G N^o : -- < x < -, 
n n 

and we will write x ~ 0. Similarly, we write y^ziiy — z^O. 
Clearly, such an x G *]Rf will be infinitesimal if and only if in the field 
of fractions *Mfi.ac, the element x~^ is infinitej^ i.e. it doesn't satisfy a 
bound of the form Ixl < n for some n G N. 



^Moreover, examining the conditions defining the notion of ultrafilter, one can 
guess that the notion of a dominant set is not intuitively so clear. In point of fact, 
the technically desirable conditions about the closure with respect to intersection 
and complement can lead to counter intuitive consequences. We would have that 
even numbers P2 or odd numbers will be dominant (but not both). Let us suppose, 
e.g., the first case and continue: even numbers in P2, i.e. the set P4 of multiples of 4, 
or its complement N\ P4 will be dominant. In the latter case, also P2 D (N\ P4), i.e. 
numbers of the form 2(2n + 1), will be dominant. In any case we would be able to 
find always a dominant set which has "1 /2 of the elements of the previous dominant 
set". Continuing in this way, we can obtain a dominant set, which is intuitively very 
"thin" with respect to its complement. To understand this idea a little better, let 
us consider that everything we said up to now can be generalized if instead of 
sequences m : N ^ R we take functions u : [0, 1] — >■ IR. In other words, instead of 
taking our indices as integer numbers, we take real numbers in [0, 1]. Then, we can 
repeat the previous reasoning considering, at each step k, subintervals of length 
2~''. Therefore, for every e > 0, we are always able to find in an ultrafilter on [0, 1] 
a dominant set A whose uniform probability P{A) < e, whereas P{[0, > 1 — e, 
even though this complement is not dominant. See [5] for a formalization of this 
idea using the notion of density of subsets of N. 

''Note that without additional assumptions, *Rfrac may contain additional more 
infinitesimals not found in *Rf; see Remark [7] below. 



8 PAOLO GIORDANO AND MIKHAIL G. KATZ 

Do infinitesimals in *]Rf correspond to ordinary null sequences? 
Theorem 5. Let [u] G *]R/, then we have that 

[u] is infinitesimal 

if and only if 



limM„ = 

n 



Proof. Let us assume that [u] is infinitesimal, then for each n G N^o, 
the set 

A„ := J G N : <Uk< - 

I n n 

is dominant. Therefore, it is infinite and we can always find an increas- 
ing sequence /c : N — N such that kn G An and kn+i > kn- For such a 
sequence we have 

Vn G N^o ■■ < < -• 

n n 

Since u G C is a Cauchy sequence, we obtain 

lim Un = lim Uk„ = 0. 

n— >+oo n— >+oo 

To prove the converse implication, we can consider that 

Vn G N^o 3N : \/k e N>Ar : - < < -. 

~ n n 

Since every cofinite set N>7v is dominant, this proves that [u] is infini- 
tesimal. □ 

As a corollary, we have that [u] ~ [v\ if and only if lim„M„ = lim„f„. 
This allows us to define the standard part mentioned above. 

Definition 6. Let [u] G *Mf, then the real number 

st([n]) := lim u„ G M (4.1) 

n— >+oo 

is called the standard part of [u] . 

Note that we have x ~ st(x) for every x G *]Rf. 

Is our extension *Mf of M still Cauchy complete with respect to some 
kind of metric extending the usual Euclidean metric on M? It is not 
hard to prove that 

:= |st(x) - st(?/)| G M Vx,?/G*Mf 

defines a pseudo-metric having the desired properties. Note that *]Rf is 
not Dedekind complete, since the set of all the infinitesimals is bounded 
but does not admit a least upper bound. 
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Remark 7. In the quotient field *Mfracj the assertion of Theorem [5] is 
not generally true, unless one considers a particular type of ultrafilter, 
called a P-point. While the existence of a free ultrafilter can be proved 
using Zorn's lemma, which is equivalent to the axiom of choice, the 
existence of a P-point cannot be proved in ZFC, that is using the usual 
axioms of set theory plus the axiom of choice. Assuming the contin- 
uum hypothesis or Martin's axiom and using transfinite induction, it 
is possible to prove the existence of a P-point. See |5] and references 
therein for more details about this foundational wrinkle^ 

5. Leibniz's law of continuity in *Rp 

To convey the full power of the idea (12. ip . we have to go back to 
Leibniz. Leibniz introduced infinitesimal and infinite quantities, and 
developed a heuristic principle called the "law of continuity", which 
had roots in the work of earlier scholars such as Nicholas of Cusa and 
Johannes Kepler. It is the principle that: 

What succeeds for the finite numbers succeeds also for 
the infinite numbers 

(see Knobloch p. 67], Robinson [321 P- 266], and Laugwitz [27]). 

Kepler had already used it to calculate the area of the circle by rep- 
resenting the latter as an infinite-sided polygon with infinitesimal sides, 
and summing the areas of infinitely many triangles with infinitesimal 
bases. Leibniz used the law to extend concepts such as arithmetic oper- 
ations, from ordinary numbers to infinitesimals, laying the groundwork 
for infinitesimal calculus. 

Of course, a modern mathematical version of this heuristic law de- 
pends on our formalization of the first word 'what' in the law of conti- 
nuity as stated above. We have already seen that this is almost trivial 
if what we mean by the Leibnizian 'what' is "continuous equalities 
between real numbers". In fact we can extend arbitrary continuous 
functions as follows. Recall that C is the space of Cauchy sequences of 
real numbers. 



The previous Theorem [S] can be easily extended to *]Rfiac if we consider fractions 
■TT for which the hmit hm„^-i_no — exists finite. It results that rr is infinitesimal 
in *Mfrac if and only the limit of this fraction is zero. This permits to extend the 
definition of the standard part function to all the fractions which are of the form ^ 
but whose limit exists finite. Finally, because of the previous Remark [7] and of the 
isomorphism *Kfrac — *'^, we can always find an infinitesimal -pj € *IRfiac which is 

not generated by an infinitesimal sequence ^-^^ ; of course, this fraction is of the 

form ^ but the corresponding ratio of sequences doesn't converge. 
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Definition 8. Let / : R'^ — t- M be a continuous function. Tlien we liave 
/ o (n^, . . . , u'^) G C for every d-tuple of Caucliy sequences u^, . . . ,u'^ ^ 
C, and we can define tlie extension */ by setting 

VK],...,[nlG*Mf. 

This gives a true extension of /, i.e. */(ri, . . . , r^) = /(ri, . . . , r^^) 
for every ri,...,rrf G M (identified with the corresponding constant 
sequences) . 

Theorem 9. Let f , g : M.'^ ^ be continuous functions, then the 
equality 

Vxi, . . . G M : f{xi,...,Xd)=g{xi,...,Xd) (5.1) 

is satisfied if and only if 

Vai,...,ad G *M/: *f{ai,...,ad) = *g{ai,...,ad) (5.2) 

Analogously, we can formulate the transfer of inequalities of the form 
/(xi, ...,Xd)< g{xu . . .,Xd)- 

Proof. The equality (15. ip implies 

{neN\fial...,ai)=gial...,ai)}=N, 

where [a^] = a^- The whole set N is dominant, and therefore (15. 2p 
follows. The converse implication follows from the fact that */ and *g 
extend / and g and from the embedding R C *Mf. □ 

Remark 10. The reader would have surely noted that some of the limi- 
tations we have presented can be avoided by generalizing our construc- 
tion further. For example, the ring *Mf is only an integral domain and 
not a field, because it is not closed with respect to pointwise inverse, be- 
cause the latter are not Cauchy sequences. Similarly, we cannot extend 
a general function / : M*^ — )■ M but only continuous functions because 
we need to ensure that the image sequence is Cauchy. However, all the 
ideas we have introduced up to now work if we replace C by the whole 
of M^. In this way, we obtain a field *R and arbitrary functions can be 
extended. In the setting of the ring *Mf only continuous functions can 
be extended and hence a continuity hypothesis has to be assumed if 
one wants to use its infinitesimals. For more details, see Goldblatt |14j . 
In the present article, we adhere to the framework of defining a new 
notion of completeness so as to add new infinitesimal points to M. We 
will motivate our definition by developing some powerful key properties 
of an infinitesimal-enriched extension of M. 



7(K],...,[n1):= {f{ul 



neN 
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Can Leibniz's law of continuity be proved for more general properties, 
e.g. for order relations or disjunctions of equality and inequality or even 
more general relations? To solve this problem, we start, once again, 
from a historical consideration. 

Cauchy used infinitesimals to define continuity as follows: a func- 
tion / is continuous between two bounds if for all x between those 
bounds, the difference f{x + h) — f{x) will be infinitesimal whenever h 
is infinitesimal, see [1]. 

Such a definition tends to bewilder a modern reader, used to thinking 
of / as being defined for real values of the variable x, but now we can 
think of f{x + h) as corresponding to *f{x + h). The function / is not 
necessarily defined on all of M, so that an extension of the real domain D 
of the function is implicit in Cauchy's construction. Therefore, we will 
start by defining such an extension of C M. We will first define the 
symbol "G„", and then define *Di in terms of En- 

Definition 11. Let n G C be a Cauchy sequence and C M, then 
(1) Un Gn D <^=^ {ri G N I u„ G D} is dominant 



Let us note that the variable n is mute in the notation m„ G^ D. 

Using this notation, our questions concerning Leibniz's law of continu- 
ity can be formulated as preservation properties of the operator *(— )f. 
In fact, as in Theorem [9l where equalities between continuous func- 
tions are preserved, we can ask whether *(— )f preserves intersections 
(i.e. "and"), unions (i.e. "or"), set-theoretic difference (i.e. "not"), in- 
clusions (i.e. "if... then..."), etc. To this end, it is interesting to note 
that a minimal set of extension properties necessary implies ultrafilter 
conditions. 

We will use a circle superscript ° in place of a star to indicate a 
general extension. 

Theorem 12. Assume that °(— ) : V{M.) — )■ P(°]R) preserves unions, 
intersections and complements, i,.e. for every A, B CW, we have 



(2) *Dr.= {[u]e*Ri\u„ G 



n 



D} 



o 



o 



o 



(AUB) 

(AnB) 

{A\B) 



o 



o 



o 



'An°B 
AX^B. 



AU°B 



Finally, let e G °M. Then 

7^e := {X C R|e G °X} 
and if e E °N, then 

K-= {X nN\x e Tie} 



is an ultrafilter on R, 



is an ultrafilter on N. 
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Proof. We need first to prove that °(— ) preserves also the empty set 
and inclusions. Indeed, °0 = °(0 \ 0) = °0 \ °0 = 0. Assume AC B,so 
that A = AnB and°A = °Ar]°B and thus °A C °B. 

lfX,Ye 7^e, then e G °X n °F = °(X n F), and hence X n F € 7^e• 
If X e We and M D F D X, then e e °X C °Y and hence Y G 7^e• 
If X C M, then R = X U (M \ X); but e G °M = °X U (°M \ °X) and 
therefore X G T^e or R \ X G T^e, and this finally proves that TZ^ is an 
ultrafilter on R because every X G TZe is not empty since °(— ) preserves 
the empty set. 

The proof that A/'e is closed with respect to intersection is direct. 
Consider N D 5 D X n N with X G 7^e; then Y := {S \ X) U X D X 
and hence Y G T^e- Therefore, F fl N = 5 because X fl N C S, and 
hence S G AC. Finally, if 5* C N, then either 5* G TZe, and thus 5 = 
5nN G A/'e, or R\5 G 7^e. In the second case, (R\5) HN = N\5 G AC. 
Up to now, we didn't need the further hypothesis e G °N. However, in 
this case, if X G 7^e, then e G °X n °N = °(X n N) 7^ and hence also 

xnN7^0. □ 

Taking, e.g., e = 1 G *Nf, yields an ultrafilter (a so-called principal 
ultrafilter, see the next Corollary [T3|) . 

The meaning of this theorem is the following: if one doesn't like the 
idea fl2.ip but wants to obtain something corresponding to Leibniz's 
law of continuity, one must face the problem that the corresponding 
extension operator °(— ) cannot preserves "and", "or" and "not" of ar- 
bitrary subsets. In Section 0, where we will introduce another idea to 
refine Cauchy's equivalence relation without using ultrafilters, we will 
see that a corresponding law of continuity holds, but only for open 
subsets, so that we are forced to define a set-theoretical difference with 
values in open sets 



where int(— ) is the interior operator. Note that the use of open sets 
and this "not" operator correspond to the semantics of intuitionistic 
logic. 

For the sake of completeness, we also add the following results, which 
represents particular cases of the previous Theorem [T2l 

Corollary 13. In the hypotheses of TheoremlT^ if 



then we have that e G R and only if IZe is the principal ultrafilter 
generated by e, i.e. 



A\B := mt{A\B) 



X c °x , (°X \ X) n R = VX C R, 



(5.3) 



7^e = {X C R I e G X} . 



(5.4) 
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Proof. Let us assume that e G M and prove the equahty ( 15.41) . If 
e e X C M, then e G °X because X C °X by hypotheses, and therefore 
X G Tie- Vice versa if e G °X, then °X = X U i°X\X) and hence 
e e X because, by hypotheses, {°X \ X) n M = and e G M. 

Finally, the converse implication follows directly from the equality 
(15^ and from Re Tie. □ 

Therefore, if the extension operation °X really extends X (first con- 
dition of (15. 3p ) adding only new non real points (second condition of 
(15. 3p ). then taking e G M we get a trivial ultrafilter. However, in our 
construction we started from a free ultrafilter; this is the case consid- 
ered in the following corollary. 

Corollary 14. In the hypotheses of Corollary [73|, let us assume that 
(°R, <) is an ordered set extending the usual order relation on the reals. 
Suppose that e G °M \ M infinite with respect to (°]R, <), i.e. 

ViV G N : e > iV 

and also that 

VATgN: e>N e G°[iV, +oo), 

then the ultrafilter Me is free. 

For example, the field *Mfrac satisfies the hypotheses of this corollary 
if we take e = tttyT- 

Proof. By our hypothesis, every interval [A^, +oo) = {x G M | x > A^} 
is in Tie, therefore [A^, +oo) fl N G Me- If X C N is cofinite, then 
N \ X C [0, X) for some iV G N and hence X D [iV, +oo) n N. From 
Theorem [T^l we have that Me is an ultrafilter, so that it is closed with 
respect to supersets, and hence X G Me- □ 

Our operator *(— )f has the following preservation properties of propo- 
sitional logic operators. 

Theorem 15. Let A, 5 C then the following preservation properties 
hold 

(1) *(AU5)f = *AfU*5f 

(2) *(An5)f = *y4fn*fif 

(3) *{A\B), = *A,\*B, 

(4) A C if and only if *Af C *5f 

(5) *0f = 

(6) *Ai = *B{ if and only if A = 5. 
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Proof. For example, we will prove the preservation of unions, the other 
proofs being similar. Take [u] E *{AU B)^, then {n \ Un E A U B} is 
dominant. If {n | m„ G A} is dominant, then [u] G *Af, vice versa, 
{n\un ^ A} is dominant and therefore it is also the intersection 

{n\un e AU B} n {n\un ^ A} = {n\un e B} , 

so that [u] G *Bf. Vice versa, if e.g. [u] G *Af, then {n | m„ G A} is 
dominant, and hence also the superset {n | m„ G A U B} is dominant, 
i.e. [u] G *(AU5)f. □ 

Example 16. Let A, 5, C C R and write e.g. A{x) to mean x G A. 
We want to see that our previous Theorem [TS] implies that Leibniz's 
law of continuity applies to complicated formulas like 

Vx G M : A{x) [B{x) and ^ . (5.5) 

In other words, we will show how to apply the previous theorem to 
show that fl5.5p holds if and only if the following formula holds 

Vx G *Mf : *Ai{x) [*Bf{x) and (*Cf(x) *Df{x))] , (5.6) 

where e.g. *Ai{x) means x G *Ai. In fact, if we assume ( 15. 5p . this 
implies that A C B and hence, by Theorem [151 *A{ C *B{. Therefore, 
if we assume *Ai{x), for x G *]Rf, from this we immediately obtain 
*B{{x). The hypotheses (15. 5p also implies that Ad C D, so that if 
we further assume *Cf(x) we also obtain that *Df{x) holds, and this 
concludes the proof of (15. 6p . Analogously we can prove the opposite 
implication. 

Remark 17. Of course, the previous example can be generalized to every 
logical formula, proceeding by induction on the length of the formula, 
but this requires the usual (simple) background of (elementary) formal 
logic. As it is well known (see e.g. |30[ [T5| [21 123]), our example further 
shows that this "more advanced" use of nonstandard analysis can be 
left only to selected readers. 

Now, the next problem is natural: what about the preservation of 
existential and universal quantifier? We have already considered the 
case of logical connectives like "and", "or", "not" without stressing too 
much the need to have a background in formal logic. This permits 
to simplify our presentation and opens this type of setting to a more 
general audience, including physicists and engineers. We wish to retain 
the same attitude also toward quantifiers. For this goal we consider two 
sets X, Y CM. and the projection px '■ X x Y ^ X , px{x, y) = x, and 
C C X X Y , i.e. a relation of the form C{x, y) with x G X and y &Y. 
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We have 



Px{C) 



{x e X \ 3z e C : X = px{z)} = 
{xeX\3yeY: C{x,y)}; 
{xeX\^i3yeY: ix,y) ^ C)} 
{xeX\yyeY: C{x,y)}. 



X\px [{XxY)\C] 



Therefore, now our aim is to prove that *(— )f preserves px{C), which 
corresponds to the existential quantifier (preservation of universal quan- 
tifier follows from this and from the preservation of difference). Only 
here we notice that, exactly as we proceeded for functions considering 
only the continuous ones, we need an analogous condition for a relation: 
what is a continuous relation C C X x Y7 To find the corresponding 
definition, we start from the idea that if / : M — t- M is continuous, then 
we expect that the relation {{x,y) E X x Y \ y = f{x)} is continuous. 
We can therefore note that the peculiarity of the definition of the ex- 
tension */ (see Definition [H]) is that continuity permits us to define */ 
on all of *Mf. Otherwise, we would always have the possibility to define 
*/ on the smaller domain 



For this reason, we start by introducing the following definition. 
Definition 18. Let X, F C M and C C X x F, then 



where dom(C) = {x E X\3y eY : C{x,y)} is the domain of C , and 
cod(C) = {y eY \ 3x E X : C{x,y)} is the codomain of C . 

Proof. We prove, e.g., the relation about the domains. If [u] G dom(*Cf), 
then there exists v such that ([m], [v]) G *C{, i.e. m„ Gn dom(C), and 



Therefore, it is the opposite inclusion that represents our idea of a 
continuous relation. 

Definition 20. In the previous hypothesis, we say that: 

(1) C is continuous in the domain iff dom(*Cf) ^ *[dom(C)]f. 

(2) C is continuous in the codomain iff cod(*Cf) 3 *[cod(C)]£. 



{[u] e*Rf\foueC}. 




dom(*Cf) C *[dom(C)]f 
cod(*a) C *[cod(C)]f, 



this means that [u] G *[dom(C)] 



□ 
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For example, in the case C = graph(/), the continuity in the domain 
says that */ is defined on the whole *Xf. Analogously, we can define 
the continuity of an n-ary relation with respect to its A;-th slot. 

Theorem 21. If X , Y C W, and f : X ^ Y, then f is continuous if 
and only if graph{f ) is continuous in the domain. 

The proof of this theorem can be directly deduced from the following 
consideration. The continuity of C in the domain can be written as 

VugC: n„G„dom(C) ^ By e *Yr. *Q{[u],y) . (5.7) 

We can write this condition in a more meaningful way if we use the 
following notation for an arbitrary property V{n): 

[V"^?7, : 'P{n)~^ <^=^ {n G N I 'P{n)} is dominant. 

For example, u„ Gn D can now be written as W^n : Un E D. Therefore, 
(15. 7p can be written as 

VugC: (V^n3y gF: C(n„,y)) By e *Yr. *Q{[u],y). (5.8) 

This can be meaningfully interpreted in the following way: if we are 
able to solve the equation 

C(m„, yn) = true 

finding a solution ?/.„ G F for a dominant set of indices n, then we are 
also able to solve the equation 

*C{{[u],y) = true 

for a solution y G *Yf. 

Using this formulation, it is not hard to prove that all the relations 
=, < and < are continuous both in the domain and in the codomain. 
An expected example of non continuous relation is x ■ y = 1 (take, 
e.g., Un '■= in (15. 8p ). This corresponds to the non applicability of 
Leibniz's law of continuity to the field property 

\/x eR: Xy^O 3y eR: x ■ y = 1, 

which cannot be transferred to our *]Rf, which is only a ring and not a 
field. 

Now, we can formulate the preservation of quantifiers: 

Theorem 22. Let X, F C M and C ^ X xY be a relation continuous 
in the domain, then 

*[px{C)], = p.x,CQ). 

That is 

*{xeX\3yeY: C{x,y)}, = {x e *Xi\3y e *Yr. *a{x,y)} . 
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As a consequence we also have 

*{xeX\\/yeY: C{x,y)}, = {x e *Xi\^y e *Yr. *Q{x,y)}. 

Proof. If [u] E *[px{C% then Un PxiC), i.e. 

V^n : M„ e X , 3y eY : C{un,y), 

that is the set of n G N satisfying this relation is dominant. This 
imphes that m„ G„ X and hence [u] G *Xf and m„ G„ dom(C), i.e. 
[u] G *[dom(C)]£. Our relation C is continuous, so that [u] G dom(*Cf), 
i.e. 

3/3G*Ff: *Cf(M,/3), 
which can also be written as 

To prove the opposite inclusion it suffices to reverse this deduction and 
use Theorem [T9] instead of the Definition 120) of continuous relation. □ 

Example 23. Let us apply our transfer theorems to a sentence of the 
form 

yaeABbeB: C{a,b) (5.9) 

showing that it is equivalent to 

Va G * 36 G *5f : *Cf (a, 6) , (5.10) 

where C C A x i? is a binary continuous relation. Assume (15. 9p and 
a G *Af. From (15. 9p we have that 

AC {a eR\3b e B : C(a,6)}, 

Therefore, by Theorem [TSl and Theorem [22l we have 

*Ai C {a G *Mf I 36 G *B{ : *C{{a, 6)} 

and we obtain the existence of a 6 G *-Bf such that *Cf(a, 6). To prove 
the opposite implication, it suffices to reverse this deduction. 

The ring *]Rf and the field *Mfi.ac can be used to reformulate profi- 
ciently several parts of the calculus. For example, a continuous function 
/ : M — )■ M is differentiable at a; G M if there exists m G M such that for 
every non zero infinitesimal h G *Mf 

3o- G *]Rf : f{x + h) = f{x) + h-m + h-a , a ^ 0, 

that is if f{x + h) is equal to the tangent line y = f{x) + /i ■ m up to an 
infinitesimal of order greater than /i, i.e. of the form h ■ a, with a ~ 0. 
Because *]Rf is an integral domain, taking a non zero infinitesimal h, 



18 



PAOLO GIORDANO AND MIKHAIL G. KATZ 



we can easily prove that such m G M is unique. Working in *Mfrac, we 
have that such m is given by 

Of course, this real number will be denoted by f'{x), so that we have 
that it is infinitely close (or, in Fermat's terminology, adequal; see [35| 
p. 28]) to the corresponding infinitesimal ratio 

_ f{x + h)- f{x) 

Reformulation is only the most trivial possibility offered by a contin- 
uum with infinitesimals, because our geometrical and physical intuition 
is now strongly supported by a corresponding rigorous mathematical 
formalism. 



6. A POSSIBLE APPROACH WITH NILPOTENT INFINITESIMALS 

There is another approach of refining Cantor equivalence relation on 
real Cauchy sequences. This approach avoiding ultrafilters. The idea 
is to compare two sequences v E C with a basic infinitesimal, e.g. 
(-) . We therefore set by definition 

u V <^=^ lim n ■ (u„ — f„) = 0. (6.1) 

n— >+oo 

In other words, using Landau's little-oh notation, the two Cauchy se- 
quences are to be equivalent if 

Un = Vn + (i) for H -^OO. 

As in the previous part of the article, we will denote the equivalence 
class of a sequence u simply by [u]. The relation defined in ( 16. ip is 
stronger than the usual Cauchy relation: 

M ~ f ^ 3 lim Un = lim Vn ='■ st([u]) G M. 

n— >+oo n— >-+oo 

It is also strictly stronger, because, e.g., the equivalence class 
with < p < 1, is a nonzero infinitesimal. For example, the infinitesi- 
mal [(^)^] is not zero, but we can think of it as being so small that its 
square is zero: = [0]. With respect to pointwise operations, we 

thus obtain a ring rather than a field. A ring with nilpotent elements 
may seem unwieldy; however, this was surely not the case for geometers 
like S. Lie, E. Cartan, A. Grothendieck, or for physicists like P.A.M. 
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Dirac or A. Einstein (see, e.g., references in [8]). The latter used to 
write formulas, if i;/c ^ 1, like 

1 



containing an equality sign rather than an approximate equality sign. 
More generally, in [7] A. Einstein wrote 

/(a;,t + r) = /(x,t) + r-^(x,t) (6.2) 

justifying it with the words "since r is very small". Let us note that 
if we apply (\6.2h to the function f{x,t) = at t = 0, we obtain = 
+ r ■ = and therefore we necessarily obtain that our ring of scalars 
contains nilsquare elements. Of course, it is not easy to state that 
physicists like A. Einstein or P.A.M. Dirac were consciously working 
with this kind of scalars; indeed, their work, even if it is sometimes 
found to be lacking from the formal/syntactical point of view, it is 
always strongly supported by a strong bridge with the physical meaning 
of the relationships being discovered. 

A difficult point in working with a ring having nilpotent elements is 
the concrete management of powers of nilpotent elements, like 

'n • • • ""n ■ 

Let us note that this kind of product appears naturally in several vari- 
able Taylor formulae. Is such a product zero or not? Are we able 
to decide effectively whether it is zero starting from the properties of 
the infinitesimals hj and the exponents ijl To be able to give an af- 
firmative answer to this, and several other questions, we restrict this 
construction to a particular subclass of Cauchy sequences, as follows. 



+ 'S^ai- \' o ( —] as n +00, (6.3) 

n"-' J 



Definition 24. We say that n is a little-oh polynomial, and we write 
n G Mo if and only if we can write 

Ur,. = r ' 

i=l 

for suitable A; G N, r, ai, . . . , a/c G M, ai, . . . , G M>o- 

Therefore, [^] C C and our previous example is generated 

by a little-oh polynomial. Little-oh polynomials are closed with respect 
to pointwise ring operations, and the corresponding quotient ring 

1 . 
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is called ring of Fermat reals. The name is motivated essentially by 
two reasons: in the ring of Fermat reals, it is possible to formalize 
the informal method used by A. Fermat to find maxima and minima, 
see |11J ; all the theory of Fermat reals and Fermat extensions has been 
constructed trying always to have a strong bridge between formal prop- 
erties and informal geometrical interpretation: we think that this has 
been one of the leading methods used by A. Fermat in his work. For 
all the proofs of this section, we refer to [T0| |8]. 

Exactly as in the previous section about the hyperreals, we have that 
the ring of Fermat reals *M is still Cauchy complete with respect to the 
pseudo-metric 

d{x,y) := \st{x) - st{y)\ eM Wx,ye'R. 

Once again, the ring 'M. is not Dedekind complete. 

It is not hard to prove that all the numbers k,r,ai,ai appearing in 
(16. 3p are uniquely determined if we impose upon them the constraints 

< fli < a2 < ■ ■ ■ < flfe < 1 (6.4) 
ai^O yi = l,...,k. (6.5) 
We can therefore introduce the following notation. 

Definition 25. If x := [u] G *M and k,r,ai,ai are the unique real 
numbers appearing in (16. 3p and satisfying (16. 4p and (16. 5p . then we set 
°x := st{x) := r, °Xi := a^, u{x) := UJi{x) := ^, := k. Moreover, 
we set 



dta :-- 




Va G 



and, more simply, dt := dti. Using these notations, we can write any 
Fermat real as 

X = °x + '^°Xi - dt^^(^) (6.6) 

i=l 

where the equality sign has to be meant in *]R. The numbers °Xj are 
called the standard parts of x and the numbers uji{x) the orders of x 
(for i = 1 we will simply use the names standard part and order for °x 
and u!{x)). The unique writing (16.60 is called the decomposition of x. 

Let us note the following properties of the infinitesimals of the form 

dta. 



dta ■ dtb = dt ah 

a + b 



dtaY = dts. Vp G M>i 
dta = Va G M<i. 
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A first justification to the name "order" is given by the following 

Theorem 26. If x E 'M. and k E N>i, then = m 'M if and only 
if°x = and u!{x) < k. 

This motivates also the definition of the following ideal of infinitesi- 
mals: 

Definition 27. If a G M U {oo}, then 

Da:= {x E'Rl^x = , uj{x) < a + 1} . 

These ideals are naturally tied with the infinitesimal Taylor formula 
(i.e. without any rest because of the use of nilpotent infinitesimal 
increments), as one can guess from the property 

a EN L>a = G •M|x"+^ = 0} . 

Products of powers of nilpotent infinitesimals can be effectively decided 
using the following result 

Theorem 28. Let hi, . . . ,hn E D^o \ {0} and ii, . . . , G N, then 

(2) h\' h'^ ^0 ^ = Efc=i TO' 

This result motivates strongly our choice to restrict our construction 
to little-oh polynomials only. 

The reader can naturally ask what would happen in case of a different 
choice of the basic infinitesimal (^)^ in the Definition fl6.ll) . Really, 
any other choice of a different infinitesimal (s„)„ will conduct to an 
isomorphic ring through the isomorphism 



i=l 



i=l 



This is the only ring isomorphism preserving the basic infinitesimals 
dta and the standard part function, i.e. such that: 

f{a-dta) = a- [Cv^)J^ 
fCx) = °f{x). 

Essentially the same isomorphism applies also to the ring defined in 
|10] . where instead of sequences, the construction is based on real func- 
tions of the form u : ]R>o IR. 
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-4dt 




7dt 







Figure 7.1. Some first order infinitesimals 



7. Order relation 



It is not hard to define an intuitively meaningful order relation on 
the ring of Fermat reals 

Definition 29. Let x, ?/ G *M be Fermat reals, then we say that x < y 
iff we can find representatives [u] = x and [v] = y such that 



For all the proofs of this section, see e.g. [T2| |8]. 

It is not hard to show that this relation is well defined on 'M. and that 
the induced ordered relation is total. This is another strong motivation 
for the choice of little-oh polynomials in the construction of the ring of 
Fermat reals. The analogous of Theorem |5] is the following 

Theorem 30. Let h G 'M, then the following are equivalent 
(1) /i G Doo, i-e. °h = 0, i.e. h is an infinitesimal 
2 Vn G N>o : < h < ^ 

The following theorem permits to decide algorithmically the order 
relation between two Fermat reals, using only their decompositions 

Theorem 31. Let x,y G 'M. If°x ^ °y, then 



Otherwise, if °x = °y, then 

(1) If uj{x) > uj{y), then x > y if and only if°Xi > 0. 

(2) If u{x) =u{y), then 

°xi >°yi x> y 
°xi <°y^ =^ X <y. 

For example, < dt < dt2 < dt^, etc. This motivates why we take ^ 
in the definition of dta- in this way the greater is the order a and the 
greater is the infinitesimal. 

The ring 'M can also be represented geometrically. 



3N G NVri >N: u„ < u 
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Definition 32. If x G 'M and 6 G ]R>o, then 

gmphs{x) := | (°a; + ^ °a;, ■ t^/^'^^), t) | < t < 6^ (7.1) 

E.g. graph^(dt2) = {(v^, t)|0 < t < 6}. Note that the values of 
the function are placed in the abscissa position, so that the correct 
representation of graph^(x) is given by the figure [TTTl This inversion of 
abscissa and ordinate in the graph^(x) permits to represent this graph 
as a line tangent to the classical straight line M and hence to have 
a better graphical picture. Finally, note that if x G M is a standard 
real, then = and the graph^(x) is a vertical line passing through 
°x = X, i.e. they are "ticks on axis". 

The following theorem introduces the geometric representation of the 
ring of Fermat reals. 

Theorem 33. If S & ^>0! then the function 

X e'Rt-^ graphs (x) C 

is injective. Moreover if x, y G *M, then we can find 6 G M>o (depending 
on x and y ) such that 

X < y 

if and only if 

yp,q,t : {p,t) G graphsix) , {q,t) G graph^iy) ^ p<q (7.2) 

that is if a point {p,t) on giaph.g{x) comes before (with respect to the 
order on the x-axis) the corresponding point {q,t) on graph^(?/). 

8. Infinitesimal Taylor formula and computer 
implementation 

What kind of functions / : M — )■ M can be extended on 'M? The idea 
for the definition of extension is natural */([«]) := [f o u] so that we 
have to chose / so that: 

(1) If M is a little-oh polynomial, then also / o m is a little-oh poly- 
nomial. 

(2) If [u] = [v], then also [fou] = [fov]. 

The second condition is surely satisfied if we take / locally Lipschitz, 
but the first one holds if / is smooth. 

Definition 34. Let / : M'^ — > M be a smooth function, then 

7([wi], . . . , M) := [/ o Ml, . . . , / o u,] W[u,], [ud] G -M. 
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Therefore, the ring of Fermat reals seems potentially useful e.g. for 
smooth differential geometry (see e.g. chapter 13 of jB]) or in some 
part of physics (see e.g. [^), where one can suppose to deal only with 
smooth functions. 

In several applications, the following infinitesimal Taylor formulae 
permit to formalize perfectly the informal results frequently appearing 
in physics. 

Theorem 35. Let x G M and / : M — ?■ R a smooth function, then 



In this case we have m = f'{x), where f'{x) is the usual derivative of 
f at X. 



Theorem 36. Let x G W^, n G N>o and / : M'^ — t- M a smooth function, 
then 



Note that m = f'{x) G M in Theorem [351 the slope is a stan- 
dard real number, and that we can use this formula with standard real 
numbers x only, and not with a generic x G *M, but it is possible to 
remove these limitations (see [TT| [9| [8]). 

The definition of the ring of Fermat reals is highly constructive. 
Therefore, using object oriented programming, it is not hard to write a 
computer code corresponding to 'M. We (see also [13]) realized a first 
version of this software using Matlab R2010b. 

The constructor of a Fermat real is x=FermatReal(s ,w,r) , where 
s is the n + 1 double vector of standard parts (s(l) is the standard 
part °x) and w is the double vector of orders (w(l) is the order uj{x) 
ii X G 'M \ M, otherwise w=[] is the empty vector). The last input 
r is a logical variable and assumes value true if we want that the 
display of the number x is realized using the Matlab rats function 
for both its standard parts and orders. In this way, the number will 
be displayed using continued fraction approximations and therefore, 
in several cases, the calculations will be exact. These inputs are the 
basic methods of every Fermat real, and can be accessed using the 
subsref, and subsasgn, notations x.stdParts, x. orders, x.rats. 
The function w=orders(x) gives exactly the double vector x. orders if 
X G *]R \ M and otherwise. 



3! m G M V/i G : /(x + h) = f{x) + h ■ m. 



(8.1) 




lil<" 
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The function dt(a), where a is a double, construct the Fermat real 
dta- Because we have overloaded all the algebraic operations, like x+y, 
x*y, x-y, -X, x==y, x~=y, x<y, x<=y, x"y, we can define a Fermat real 
e.g. using an expression of the form x=2+3*dt (2) -l/3*dt (1) , which 
corresponds to x=FermatReal ( [2 3 -1/3] , [2 l],true). 

We have also realized the function y=decomposition(x) , which gives 
the decomposition of the Fermat real x, and the functions abs(x), 
log(x), exp(x), isreal(x), isinf initesimal (x) , isinvertible(x). 

The function plot(t,x) shows the curve (17. ip at the given input t 
of double. 

The ratio x/y has been implemented if y is invertible. Finally, the 
function y=ext (f , x) , corresponds to '/(x) and has been realized using 
the evaluation of the symbolic Taylor formula of the inline function f . 

Using these tools, we can easily find, e.g., that 

sinfdtg + 2dt2) , ,1 5 , 

dta + 2dt2 - -die + -dt 



cos(-dt4 -4dt) 2 5 6 

This corresponds to the following Matlab code: 
» x=dt(3)+2*dt(2) 
X = 

dt_3 + 2*dt_2 
» y=-dt(4)-4*dt(l) 

y = 

-dt_4 - 4*dt 
» g=inline ( ' cos (y) ' ) 

g = 

Inline function: g(y) = cos(y) 
>> f =inline ( 'sin(x) ' ) 
f = 

Inline function: f(x) = sin(x) 
» decomposition (ext (f ,x) /ext (g,y) ) 
ans = 

dt_3 + 2*dt_2 + l/2*dt_6/5 + 5/6*dt 

The Matlab source code is freely available under open-source licence, 
and can be requested to the authors of the present article. 



9. Leibniz's law of continuity in 'M 

Is a suitable form of the Leibniz's law of continuity provable in the 
ring of Fermat reals? The first version is the transfer for equality and 
inequality, that can be proved proceeding like in Theorem [9l 
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Theorem 37. Let f , g : Mf^ ^ be smooth functions, then it results 

Vxi, . . . , G M : /(xi, ...,Xd)= g{xi, ...,Xd) 

if and only if 

Vxi, . . . , Xd G 'M : 'f{xi, ...,Xd) = *g{xi, ...,Xd). 

Analogously, we can formulate the transfer of inequalities of the form 
/(xi, . . . ,Xd) < g{xi, . . . ,Xrf). 

Now, we can proceed as for *Mf. We firstly define the extension *U 
of a generic subset t/ C M. 

Definition 38. Define the set of little-oh polynomials Uq [^] as in 
Definition [21] but taking sequences n : N — ?■ f/ with values in U and 
such that °[u] := lim„_^+ooW„ G U. For u, v ^ Uo [^] define m ~ for 
Un = Vn + o (i) as n -> +00 and 'U := Uo [^] / ~. 

li i : U M is the inclusion map, it is easy to prove that its Fermat 
extension 'i : *U ^ *M is injective. We will always identify 'U with 
'i{*U), so we simply write 'U C 'M. According to this identification, 
if U is open in M, we can also prove that 

= {x e'R\°x eU}. (9.1) 

Because of our Theorem [T2] we must expect that our extension operator 
'(— ) doesn't preserve all the operators of propositional logic like "and", 
"or" and "not". To guess what kind of preservation properties hold for 
this operator we say that the theory of Fermat reals is strongly inspired 
by synthetic differential geometry (SDG; see, e.g., [211 ED [I])- SDG 
is the most beautiful and powerful theory of nilpotent infinitesimals 
with important applycations to differential geometry of both finite and 
infinite dimensional spaces. Its models require a certain knowledge of 
Topos theory, because a model in classical logic is not possible. Indeed, 
the internal logic of its topos models is necessarily intuitionistic. Fer- 
mat reals have several analogies with SDG even if, at the end it is a 
completely different theory. For example, in 'M the product of any two 
first order infinitesimals is always zero, whereas in SDG this is not the 
case. On the other hand, the intuitive interpretation of Fermat reals is 
stronger and there is full compatibility with classical logic. 

This background explain why we will show that our extension op- 
erator preserves intuitionistic logical operations. Even if the theory of 
Fermat reals can be freely studied in classical logicQ, the "most natural 
logic" of smooth spaces and smooth functions remains the intuitionistic 
one. We simply recall here that the intuitionistic Topos models of SDG 



'More generally, without requiring a background in formal logic. 
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show formally that L.E.J. Brouwer's idea of the impossibility to define 
a non smooth functions without using the law of excluded middle or 
the axiom of choice is correct. 

Because we need to talk of open sets both in M and in 'M we have 
to introduce the following 

Definition 39. We always think on 'M the so-called Fermat topology, 
i.e. the topology generated by subsets of the form *U C 'M. for U open 
in M. 

Theorem 40. Let A, B he open sets of M, then the following preser- 
vation properties hold 

(1) '{AUB) = 'AU'B 

(2) '{Ar\B) = Mn'B 

(3) •int(A \B) = int(M \ 'B) 

(4) A C 5 if and only if 'A C 'B 

(5) -0 = 

(6) M = 'B if and only if A = B 

Proof. We will use frequently the characterization (19. ip . To prove ([1]) 
we have that x e '{AU B) iS x e 'R and °x e AU B, i.e. iff °x G A 
or °x & B and, using again (19. ip . this happens iff x G *A or x G 'B. 
Analogously, we can prove ([2]). We firstly prove (jlj). li A ^ B and 
X G 'A, then °x E A and hence also °x E B and x G 'B. Viceversa if 
'A C *B and a E A, then °a = a so that a G 'B, that is °a = a E B. 
To prove we have that x G 'mt{A \ B) iff °x G int(74 \ 5), i.e. 
iff {°x - 5,°x - 8) C A \ 5 for some 5 G M>o. From ^ we have 
— 5, °x — (5) C 'A and x G *{°x — 5,°x — 5). Finally, a generic 
y G *(°a; — 5, °a; — 5) cannot belong to *B because, otherwise, °y G 
{°x — 5,°x — 5) r\ B which is impossible. Therefore, x is internal to 
*A \ *B with respect to the Fermat topology. The proofs of and 
are direct or follow directly from (jl]). □ 

Example 41. Using the previous theorem, we can prove the transfer 
of the analogous of ( 15. 5p . but where we need now to suppose that A, 
B, C are open subsets of M. Therefore, we have 

Vx G M : A{x) [B{x) and {C{x) =^ D{x))] 

if and only if 

Vx G 'M : M(x) ['B{x) and {'C{x) 'D{x))] . 

Once again, we don't strictly need a background of intuitionistic 
logic to understand that the preservation of quantifier for the Fermat 
extension *(— ) must be formulated in the following way 



28 



PAOLO GIORDANO AND MIKHAIL G. KATZ 



Theorem 42. Let A, B he open subsets o/M, and C be open in Ax B. 
Letp : {a,b) ^ Ax B a & A be the projection on the first component. 
Define 

•C:={(a,/3)|(°a,°/3)GC} 
3,{C) ■.= p{C) 

Wp{C) := int {A \ 3p (int {{A x B)\ C))) . 

Then 

• [3,(C)] = 3.,CC) 
•[V,(C)]=V.,rC). 

That is 

'{aeA\3beB: C(a, b)} = {a e 'A\3b e 'B : 'C{a, b)} 
'{aeA\\/beB: C{a, 6)} = {a G M | V6 G '5 : 'C{a, b)} . 

Proof. The preservation of the universal quantifier foUows from that of 
the existential quantifier and from property ([3]) of [40], so that we only 
have to prove * [p(C)] = 'p{*C). Consider that the projection is an 
open map, so that p{C) is open because C is open in Ax B. Therefore 
X e ' [p{C)] iff °x G p{C), and this holds iff we can find (a, 6) G C 
such that °x = p{a,b) = a E A. Therefore, *p{x,b) = [{p{xn,b))^] = 
[{xn)n] = X and (x, b) G 'C because b) = (a, b) G C. This proves 
that • [p(C)] C *p(*C). Vice versa, if x G *p{'C), then we can find 
(a, f3) G *C such that x = *p{a, (3) = a. Therefore, {°a, °(3) G C and 
= °a = °x. This means that °x G p(C), which is open and 
hence x G *p{C). □ 

Example 43. Using the previous theorem, we can prove the transfer 
of the analogous of Example [231 but where we need now to suppose 
that A, B are open subsets of R and C is open in A x B. Therefore, 
we have 

\/aeA3beB: C{a,b) 

if and only if 

Va e'A3be'B : 'C{a,b). 

The theory of Fermat reals can be greatly developed: any smooth 
manifold can be extended with similar infinitely closed points and the 
extension functor *(— ) has wonderful preservation properties that gen- 
eralize what we have just seen on the (intuitionistic) Leibniz's law of 
continuity in 'M. Potential useful applycations are in the differential 
geometry of spaces of functions, like the space of all the smooth func- 
tions between two manifolds. 
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10. Conclusion 

We started with the idea of refining the equivalence relation among 
real Cauchy sequences so as to obtain a new infinitesimal-enriched con- 
tinuum. We have developed this idea in two directions. The first di- 
rection takes one toward the hyperreals, and we tried to motivate the 
choices one must make to arrive at a powerful theory. On the other 
hand, we saw that the intuitive interpretation of such choices is some- 
times lacking. The second idea is intuitively clearer but surely formally 
less powerful. The two ideas serve different scopes because they deal 
with different kinds of infinitesimals: invertible and nilpotent. 
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